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Abstract 



We describe a Monte Carlo event generator for the simulation of baryon- and 
lepton-number violating processes at supercolliders. The package, HERBVI, is 
designed as a hard-process generator interfacing to the general hadronic event 
simulation program HERWIG. In view of the very high multiplicity of gauge 
bosons expected in such processes, particular attention is paid to the efficient 
generation of multiparticle phase space. The program also takes account of 
the expected colour structure of baryon-number violating vertices, which has 
important implications for the hadronization of the final state. 
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Program summary 



Title of program: HERBVI 

Program obtainable from: gibbsOdxcern . cern . ch 

Prgramming language used: FORTRAN 77 

Memory required: approx. 120 kbytes 

Number of bits per word: 32 

Subprograms used: HERWIG version 5.7 or later 

Number of lines in distributed program: 4000 

Keywords: baryon and lepton number violation, Monte Carlo techniques, high particle 
multiplicites 

Nature of physical problem: simulation of baryon and lepton number violating pro- 
cesses with an associated large gauge boson multiplicity. 

Method of solution: generation of phase space configurations using an efficient algo- 
rithm, coupled to a phenomenological model of the physical processes under study. 

Typical running time: 2-3 hours for a 10*^ event simulation (DEC ALPHA 4000), 
dependent on the process simulated 

Unusual features of the program: none. 



1 Introduction 



It has been known for some time that topologically non-trivial gauge field configurations 
induce baryon- and lepton-number violating (BLNV) vertices in the electroweak sector of 
the standard model |l|, §]. More recent investigations of these processes suggest that such 
processes may be observable at supercollider energies. An extensive review can be found 
in i. 

However, at present it is not possible to make definite predictions even about the rate 
of such processes. What is possible is to take the main features of the calculations and 
use these characteristics to form a phenomenological model of BLNV processes. The aim 
of such a model is to predict possible experimental avenues that can be explored, with 
the possibility of improved theoretical calculations in the future allowing the model to be 
refined. One such model has been presented in [^, ^, and in section 2 we review this model. 

In section 3 we describe the implementation of the model as a Monte Carlo simulation, 
HERBVI. This package simulates hard scattering processes involving BLNV, and has been 
designed to take into account the very high particle multiplicities associated with BLNV 
processes. In addition, the formation of the extra baryons or antibaryons from the BLNV 
process is dealt with. 

Finally, in section 4 we make some concluding remarks. In an appendix we describe the 
use of the HERBVI package by detailing how it interfaces with the general Monte Carlo 
simulation program HERWIG 0, |^. 

2 Simulating baryon number violation 

In this section, the model of |^ is briefiy outlined. One of the general properties that 
has to be addressed is the parametrically large multiplicity of gauge bosons produced in 
BLNV interactions p, |^. Consequently, we discuss in some detail the method of for 
the efficient Monte Carlo generation of many-body phase space configurations. 

2.1 General properties 

The general features of instanton-induced BLNV taken into account by HERBVI are: 

• Quarks and leptons of all families attached to the BLNV vertex, with total change in 
baryon and lepton number 

AB = AL = -3. (1) 

• A steeply rising cross section, the rise being cut off at some energy due to unitarity 
constraints. 

• A parametrically large boson multiplicity, ub, with a typical value being = 30 ~ 
1/aw- 
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The BLNV process is assumed to be initiated by a valence quark-quark collision. Thus the 
reaction is 

q + q^lq + ?>l + nBW {Z)+nHH. (2) 

where the two incoming quarks and one outgoing antiquark are first-generation, and three 
outgoing antiquarks of each of the second and third generations are produced. The matrix 
element is taken to be constant, so that the distributions of produced particles are deter- 
mined by phase space. Two options for the gauge boson multiplicity ns are provided in 
the program: either a fixed multiplicity, or a distribution modelled on the leading-order 
BLNV matrix element. Alternatively, a user-supplied routine can be used. The Higgs bo- 
son multiplicity uh is assumed to be tied to that of gauge bosons, as discussed in sect. 3.5 
below. 

For the energy dependence, a simple step-function of the quark-quark centre-of-mass 
collision energy is assumed: 

(v^) = a^e (VI - V^) (3) 



a 

The reasoning behind a model of this form has been discussed in detail in [^]. Again, a 
user-supplied subroutine may substituted if desired. 

The threshold energy parameter a/Iq used in modelling the cross section is expected 
to be of the order of the characteristic energy of these processes, the sphaleron mass [pT[ 



Msp = VGTim^/aui ~ 18 TeV. The appropriate value of the cross section parameter (Tq is a 
matter of great theoretical uncertainty. However, if one wishes to investigate the properties 
of events rather than their expected number, the magnitude of (Jq is not relevant. For 
reference, we note that the unitarity bound for S-wave processes is given by 

(Ju (s) = — (4) 
which is of the order of tens of picobarns at the sphaleron energy M, 



sp- 



in addition to BLNV events, HERBVI is able to generate the main expected background 
process, non-perturbative B- and L-conserving multi-W production. The reaction is as- 
sumed to be 

gi + ^2 ^ 93 + 94 + it-bW (Z) + HhH. (5) 

This process is also expected to occur at large boson multiplicities, ~ with an 

energy scale ~ mw/aw A detailed discussion, based on the argument that at high energies 
electroweak theory has a large, approximately constant total inelastic cross section, can 
be found in |T^. Clearly, these events will be like BLNV ones, but without any fermion 
number violation. The predominance of boson over fermion production means that the 
experimental signatures of the two types of events will be similar. It is also expected that 
the two processes will have similar threshold behaviour and cross sections. In HERBVI they 
are therefore modelled in the same way, apart from the difference in fermion production. 



2.2 Phase space generation 

The main requirement of a Monte Carlo event generator for simulating BLNV processes is 
to generate configurations in phase space for high particle multiplicities at high efficiency. 
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Therefore, careful consideration of the method used to generate phase space configurations 
is required. 



One approach is to generate n-body phase space by factorising the process into n — 1 
two-body decays. This is the method used in, for example, the phase-space generator 
GENBOD |13|. However, this method results in non-uniform phase space distributions, and 
in particular becomes very inefficient at high multiplicities. In the efficiency of this 
method is quoted as being less than 1% for n > 17. 



The MAMBO algorithm is designed to overcome this problem. The aim is to 
generate efficiently N momenta pf such that the total momentum 



N 



P' = Ep' (6) 

i=l 

is equal to the total energy and momentum required. In the rest frame of the interaction, 
this is simply 

po = yi, P = 0. (7) 

In order to develop an efficient algorithm, one may divide the process into two steps. 
Firstly, momenta fcf are generated, with the only constraint being the mass of each 
particle, kf = rrif. Obviously, their sum will not satisfy the constraints. The second step 
is to transform these momenta, using some combination of a Lorentz boost and scaling 
transformation, so that the resultant momenta do satisfy the constraints. Note that we 
may have different scaling transforms for the space-like and time-like components of the 
momenta. 

The key to the MAMBO algorithm is the scaling transformion used. First, assume that 
a suitable Lorentz transformation, acting on the to produce , has been performed so 
that the system is now at rest. Therefore 

N 

Q' = T.^^ (8) 

i=l 

such that Q = 0. The scaling transform is now performed by setting 

gf^xgf, x = VI/Qr (9) 

Therefore, the total energy is now but the particle masses have been rescaled to become 
xrrii. The last step is to now find a ( such that 

E(CV + -^?) (10) 

i=l 

and then make the replacement 

Pi = a 



P- = ^/eq^ + ml (11) 
The Monte Carlo weight associated with this procedure is 



i=l 
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It can be shown that x < 1 imphes w < 1. Therefore, by enforcing the condition x < 1 
and rejecting events with w greater than a uniform random number between and 1, we 
can use this algorithm to obtain an unweighted distribution. 



The test results described in |10] indicate that, for the typical configurations generated 
by HERB VI, the efficiency is in the region of 30% for ~ 30, falling to ~ 6% at the upper 
end of the range, ~ 120. 



3 The HERBVI package 

The HERBVI package is designed as an 'add-on' hard process generator for the Monte Carlo 
event generator HERWIG 0] . The latter program is a general purpose generator for high 
energy processes involving hadrons. In order to explain how the two packages relate to each 
other, we now give a quick review of the various stages of event generation in HERWIG. 

The first stage is the generation of the hard process under study. This can be either 
one of the processes already within HERWIG, or alternatively a separate one such as those 
supplied by HERBVI. The kinematics and parton fiavours are generated according to the 
matrix elements under study, within appropriate phase space limits. 

Next, the package simulates the perturbative evolution of the incoming and outgoing 
partons, by generating parton showers. One of the features of HERWIG is that it correctly 
treats QCD coherence effects, which lead to an ordering of the emission angles of the partons 
emitted in the showers. This ordering results in a decrease in the transverse momentum 
as one moves away from the hard process towards either the initial state beam particles 
or the final state products. The showering proceeds until either the transverse momentum 
of emitted partons is below a cut off (final state showers) or when the evolution scale is 
small enough to be cut off by a scale associated with the structure function (initial state 
showers) . 

The end of the showers represents the point where non-perturbative effects lead to 
hadronisation. HERWIG uses a cluster model W^ to describe this process. Gluons in the 



parton shower are split non-perturbatively to qq pairs, and colour singlet pairs are joined 
together to form colourless hadronic clusters. These clusters then decay into hadrons, a 
process simulated using a simple phase space model. The decay of unstable hadrons into 
stable products is then modelled using particle data tables based on experimental mea- 
surements. Heavy quark decays are treated by regarding them as new sources of coherent 
radiation and repeating the process. 

With this sketch of the HERWIG Monte Carlo generator in mind, we now turn to the 
description of the HERBVI package. Essentially, this models BLNV by generating the 
appropriate hard process for HERWIG, although the presence of baryon-number violation 
requires careful treatment when hadronisation is modelled. The package is described by 
considering various routines within it, in an order corresponding roughly to the order in 
which they are used by HERWIG. 



4 



3.1 Init ialisat ion 



HERBVI contains an initialisation procedure designed to make the changing of default 
values as easy as possible. Without any modification, the package calls the initialisation 
routines at the start of the generation of the first event. 

SUBROUTINE HVHINI - Initialise event generation variables 

This routine performs the initialisation of the variables used in event generation. It 
must be called if certain parameters are changed by the user. The parameters available to 
the user are, at present, the components of the CKM matrix (stored in the array HVCKM), 
the cross section parameter HVPLCS, and the limits on the quark-quark subprocess energy- 
squared s. The lower limit on s corresponds to the threshold energy of the model discussed 
in section 2. The limits for s generation are set by the variables SLMIN and SLMAX. These 
variables are the logarithms of the ratios SLOWER = sq/s and SUPPER = Smax/s, where s is 
the total centre-of-mass energy squared. 

SUBROUTINE HVINIT - Main initialisation routine 

This routine is only executed once per program run, and does nothing if called again. 
All the parameters and flags are set to default values within this routine. 

Should the user wish to change any default values at the start of the run, then a call 
should be made to HVINIT first, and the parameters then changed. A call will then be made 
automatically to HVHINI at the start of the first event to recalculate derived parameters. If 
the user wishes to change parameters during a run, another call to HVHINI must be made 
after the changes. 

3.2 Event generation 

The Monte Carlo program HERWIG uses the HERBVI package for process codes IPROC= 
7000 to 7999. We now examine in detail the various steps the package performs for event 
generation in this case. 

SUBROUTINE HVHBVI - Main interface routine 

This routine provides the main event generation interface to HERWIG. Its purpose is to 
translate HERWIG variables into the appropriate subroutine calls. This is achieved using 
the HERWIG variables IPROC and GENEV. The first of these determines the type of process 
to be generated. The values used in HERBVI version 1.0 are listed in Table |l]. 

The variable GENEV determines the type of routine to be called. In order to describe 
its use, we recall first how it functions within the HERWIG package. In HERWIG, a hard 
process generation routine is called twice. The first call, with GENEV set to .FALSE., 
requests the generation of an event weight by the hard process routine. The second call, 
with GENEV now . TRUE . , asks for event generation to be completed and the event record to 
be written into the standard HEPDATA common block. The HERBVI package is structured 
in a different manner, with separate routines for weight generation (GENEV=. FALSE.), and 
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HERBVI process codes 



Code Function 



7000-7099 BLNV process 
7100-7199 Multi-W process 
7200-7299 Random selection 



Table 1: HERBVI version 1.0 process codes. The random selection is of any code in the 
range 7000 to 7199 with equal probability. 

writing the event data (GENEV=.TRUE.). The event parameters are passed between the two 
sets of routines using common blocks. 

The other function of HVHBVI is a call to HVHINI, in order to initialise variables. This 
only occurs during the first execution of the routine in each run. 

Note that the dummy version of HVHBVI must be deleted from the HERWIG package to 
allow use of HERBVI. 



3.3 Weight generation 

There are two routines to generate event weights, to cover the BLNV and multi-W cases. 
We shall discuss the BLNV case first, and then highlight the differences for the B- and 
L-conserving routine. 

SUBROUTINE HVHEVT - BLNV weight generation 

This routine generates an event weight according to the BLNV model. The first function 
of the routine is to select the incoming partons, and generate the centre-of-mass energy of 
the hard process using HVSGEN. Note that the only incoming partons allowed at present 
are u or d quarks. This restriction is justified since only collisions of protons and/or heavy 
ions are envisaged, and the high threshold energies of BLNV processes correspond to a 
high momentum fraction x even at supercollider energies. Consequently valence quark 
distributions are the only significant source of incoming partons. 

The next distribution generated is that of the bosons, accomplished using HVRBOS. Two 
multiplicity distributions are provided by this routine: a constant or a leading-order type 
of distribution, selected by means of the logical flag HVC0NT(7). These choices correspond 
to the two cases discussed in [Q. Alternatively, it is straightforward for users to provide 
their own distributions. More details are given in sect. 3.5 below. 

The fermions produced in the hard process are now selected according to a simple phase 
space model. One member of each left-handed doublet has to be chosen. If the masses of 
the two particles are mi and m2, then the probability of choosing particle 1 instead of 2 is 




(13) 
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for massive particles, where (3 is estimated using 

/3 = | ^"° ■ (14) 

^ V s — UBniB 

and Ki is the modified Bessel function of the second kind. For the case mi = 0, Eq. (0) 
has to be replaced by 

p = (1 + m2(3K^ (^2^))"' . (15) 

A good approximation, used in HERBVI, is to set p = ^ except when selecting from Tb 
doublets. The error caused by this approximation is negligible, provided nBrdBl^T^ < ^■ 

As pointed out in 0], calculations of fermion number violation by instanton mediated 
processes deal with gauge-eigenstate down-type quarks, rather than the physical mass eigen- 
state particles. Therefore, the outgoing quark distribution has to be projected through the 
CKM matrix. This is achieved by a call to HVHPQM, provided the control flag HVC0NT(3) is 
. TRUE . . 

Now that the fermion distribution has been fixed, charge conservation for the whole 
process is enforced by a second call to HVRBOS, and the structure of the event stored for 
later use. The x distributions of the two incoming partons can now also be generated, by 
a call to HVSFUN. 

The last step, if requested by changing the value of the flag HVCONT(l) to .FALSE., is 
to generate the MAMBO weight for the phase space distribution. The anticipated normal 
use of the package is for an unweighted distribution to be generated during the event data 
phase of operation, but this can be overridden should the user wish to use a weighted 
distribution. The main use we foresee for this option is the generation of distributions 
corresponding to non-trivial matrix elements. 

SUBROUTINE HVHMWE - Multi-W weight generation 

The only fundamental difference between this routine and HVHEVT is in the nature and 
number of the outgoing fermions, since B and L are conserved in this case. Again, only u 
and d quarks are allowed as the incoming partons, qi and q2 in Eq. (^), and here the same 
applies to the outgoing partons, and q^. 

3.4 Event data 

The event data insertion routines are executed if the event weight has been accepted by 
HERWIG. As for the weight generation phase, there are separate routines to cover the two 
types of processes. Again, we describe the BLNV case first and then focus on the differences 
in the multi-W case. 

SUBROUTINE HVHGEN - Generate BLNV event data 

This routine writes the event data into the HEPDATA common block, and sets HERWIG 
variables as appropriate. If the flag HVCONT(l) is .TRUE., its default value, then the 
MAMBO algorithm is used to generate an unweighted phase space configuration. 
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The HERWIG book-keeping is performed by a call to HWEONE. which sets up a 2 ^ 1 
subprocess. The internal HERWIG pointers are then modified to take account of the final 
state particles. In particular, the colour connection pointers for the fermions are set in such 
a way that the outgoing light antiquark is connected to the incoming quarks, and the three 
outgoing quarks of each higher generation are connected to each other, so as to form three 
distinct colour singlets. 

SUBROUTINE HVHMWG - Generate multi-W event data 

The main difference between this and HVHGEN is the colour structure of the final state. 
The colour connections to be made between the partons are now those appropriate to 
quark-quark elastic scattering. For this purpose, the routine HVETWO is used to set up a 
2 — 2 + n process. This is a modification of the HERWIG routine HWETWO, which generates 
a 2 ^ 2 process with the correct colour pointers. 



3.5 Generation of distributions 

There are a number of routines within HERBVI to generate the various distributions. They 
are listed here together for reference. Note that these routines obtain random numbers by 
calls to random number generation routines within the HERWIG package (named HWR***). 

SUBROUTINE HVHPQM - Project quarks through CKM matrix 

The elements of the CKM quark mixing matrix are used to project the down-type quarks 
from their gauge eigenstates onto mass eigenstates. Note that if for any reason the CKM 
parameters are changed by the user, the initialisation routine HVHINI must be called to 
reset the parameters used by this routine. 

FUNCTION HVRBIN - Binomial distribution 

Generation of the binomial distribution 

The distribution is generated by first generating two uniformly random variables u and v, 
such that v? <V {v/u), and then setting n =l}i7{v/u). 

The main use of this routine is in the generation of the 7/2^° distribution by splitting 
the number of neutral bosons N into photons and Z^^s with probabilities p — sin^ 9yj and 
1 — p = cos^ 9u, respectively. 

SUBROUTINE HVRBOS - Boson number distribution 

This routine has two purposes, determined by a flag when it is called. The first is to 
generate the total number of gauge and Higgs bosons in the final state. The second function 
is to ensure charge conservation by adjusting the relative numbers of charged bosons. There 
are three control fiags HVC0NT(7,8,9) for this routine. 

The flag HVC0NT(7) controls the number of gauge bosons, ub- If .TRUE., the default. 
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then ub = 30. Otherwise, the number is selected according to a parametrisation of the 
multiphcity distribution. The defauh parametrisation is of the form 



s 



V{n,s) = expi-^—^\ (17) 

where 

a = -21.497D0 + 0.590D-2 X 

(or< 

b = 0.483D-1 + -0.854D-6 X (18) 

are parameters found by fitting to the resuhs of phase space calculations using the leading 
order matrix element for instanton-induced BLNV |16 . 



GeV 



The routine can be replaced by a user-supplied version, HURBOS, by setting HVC0NT(8) 
to .FALSE.. Note that both functions of the routine have to be catered for, and that the 
dummy version of HURBOS must then be deleted from the HERBVI package. 

The final fiag HVC0NT(9) controls the number of Higgs bosons generated. The default is 
to generate according to a simple model, namely that the probability of generating a Higgs 
boson is 1/16 of the probability of emitting a gauge boson from the BLNV vertex. This 
gives a reasonable agreement with the phase space calculations performed in the energy 
region of interest with mn = 300 GeV. Note that the leading-order instanton calculation 
yields 

which, for the region of interest E ~ M^p, also predicts the ratio uh/ub ~ 1/16. The 
phenomenology of Higgs boson production in BLNV interactions has not yet been considered 
in detail, essentially because of this expected strong suppression relative to gauge boson 
production. Setting HVC0NT(9) = .FALSE. turns off Higgs boson generation completely. 

Once the charge distribution of the particles has been fixed, the numbers of 7 and 
bosons are determined using HVRBIN as outlined above. 

SUBROUTINE HVSFUN - Parton momentum fraction generation 

This routine generates the momentum fractions x of the incoming partons. It requires 
the centre-of-mass energy of the hard subprocess to be already determined. A weight is 
generated according to the parton distribution functions of the colliding hadrons. 

SUBROUTINE HVSGEN - Hard process energy 

Generation of the total centre-of-mass energy s for the hard subprocess. The default is 
to generate s according to the distribution ds/s. Setting the fiag HVC0NT(4) to .FALSE. 
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changes this distribution to ds/s^. A user-supphed routine, HUSGEN, can be called instead 
of HVSGEN if desired by setting HVC0NT(6) to .FALSE.. 



Note that generation of a ds/s distribution is equivalent to postulating a flat parton- 
level cross section. This can be seen by noting that the convolution of the parton level 
cross section with the proton structure functions, 

pi pi PS 

(7 (s) = / dxi / dx2 / ds/i (xi) /2 (0:2) S{s - X1X2S) a (s) , (21) 
Jo Jo Jo 

gives 

.w = /;^:.m/;|.awI/.(i), (22) 

where r — s/s. Specialising to the case of a threshold parton- level cross section 

a (s) = aoe (y§ - ys^) , (23) 

we find 

where now tq — sq/s. Therefore, to generate a flat parton- level cross section, we generate 
r according to dr/r (corresponding to ds/s) and x according to dx/x, then reweight using 
the structure functions zfi{z) at z — x and z — rjx. 

Recall that the range of s is limited by the parameters SLMIN and SLMAX, which are the 
logarithms of SLOWER and SUPPER where 

SLOWER < - < SUPPER . (25) 
s 

The default values correspond to the range 

0.24 < - < 0.8 (26) 

s 

which gives between approximately 50% and 90% of the total energy. The upper limit 
is somewhat arbitrary; the rapid decrease of the structure functions at large x means that 
the contribution from the upper end of the energy spectrum is negligible. This limit has 
been introduced to increase the efficiency of the Monte Carlo. 

The user can change the variables SLMIN and SLMAX to alter the range of s generation 
at any time, although changing these parameters in the middle of simulation requires a call 
to HVHINI to enable the derived parameters to be recalculated. 



3.6 Hadron formation 



The colour structure of BLNV processes requires special treatment during the later stages 
of event processing by HERWIG. As explained above, there are three colour-singlet three- 
quark (antiquark) vertices associated with the process, the first of which comprises the 
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two incoming quarks and one outgoing antiquark. The other two correspond to outgoing 
antiquark triplets. 

These vertices do not occur in any baryon-number conserving process, and therefore 
must be treated by the HERBVI package. However, the set of three partons forming a 
singlet in the hard subprocess will not necessarily be those forming the 'extra' baryon as a 
result of BLNV, because the first gluon emitted by any of these partons will carry with it 
the colour connection to the BLNV vertex. Therefore, the parton cascade within HERWIG 
must be allowed to occur before the colour connections are made. The point in HERWIG 
where this occurs is at the start of the clustering phase of the simulation, in the routine 
HWCFOR. 

SUBROUTINE HVCBVI - Find unpaired partons after BLNV 

This routine is called at the start of the HERWIG routine HWCFOR, if a BLNV event has 
been generated. It searches through the event record to locate all partons triplets connected 
with the BLNV vertices. The end product of the routine is to produce, for each vertex, 
a quark and associated diquark (or antiquark and antidiquark, as appropriate) that are 
colour connected. 

HERWIG is at present unable to incorporate heavy quarks into diquarks. This means 
that a vertex such as bbu requires special treatment. The approach used is to form a cluster 
from two of the BLNV particles, one of which is a heavy antiquark Q. The cluster is then 
split, producing an extra qq pair. A mesonic colour singlet is formed from qQ leaving the 
hght antiquark q colour-connected to the BLNV vertex. This approach fails when there is 
insufficient energy to form the extra qq pair. In this case the event is rejected by calling 
the HERWIG error handler HWWARN with a negative error code, allowing the event to be 
rejected in a controlled manner. Typically, about 1% of events have to be rejected in this 
way. 

At the end of this routine, all the BLNV colour connections have been correctly made 
and the HERWIG package can now complete the process of event generation. Note that, 
in order to generate BLNV events, the dummy version of HVCBVI must be deleted from 
HERWIG. 



3.7 Analysis of events 

The HERBVI package contains a routine to aid in the analysis of BLNV events. Whilst this 
code is not a necessary part of the package, it does aid in the identification of simulated 
particles direct from the BLNV vertex as opposed to those formed as decay products of 
other unstable particles. 

SUBROUTINE HVANAL - Analyse final state products 

The event record is scanned to identify the particles from the BLNV vertex. Particles 
are fiagged as 'primary' from the BLNV vertex if they are either (anti)leptons from the 
vertex itself, photons from the vertex, or (anti)baryons formed from the colour triplet at 
the vertex. In addition, some other useful parameters such as the rapidity, pseudo-rapidity 
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HVADAT contents 


i\L 1 dy 




IDHNUM 

PSRFSP 

RAPFSP 

PTFSP 

IDHFSP 

ISFSP 


Pointer to entry in HEPEVT 
Pseudo-rapidity of particle 
Rapidity of particle 
Transverse momentum of particle 
Type of particle (PDG code) 
TRUE for particle from BLNV vertex 



Table 2: Contents of the HVADAT common block. Each array listed has one entry for each 
particle in the final state after the HERWIG simulation phase. 

and transverse momentum of each particle are also calculated. The data are stored in the 
common block COMMON HVADAT. 

The contents of HVADAT are complementary to HERWIG output contained in the stan- 
dard HEPEVT common block, and are listed in Table |. Note that the file HERBVIIO.INC 
contains the variable declarations for HVADAT. The variable FSPPTR, which is also stored in 
this common block, contains the total number of final state particles found and processed 
by HVANAL. 

3.8 Program summary 

For convenience we list here the main routines, control flags and default settings. The rou- 
tines are listed in Table ^ Other auxiliary routines, not listed here, all have names starting 
with either HV**** or BVP***. This convention has been followed to avoid clashing with 
other software. These routines are associated with the internal workings of the HERBVI 
package and do not concern us here. 

The control flags are listed in Table ^. The default value in all cases, as set by HVINIT, 
is . TRUE . . The default values of physical parameters are in general taken from the cor- 
responding HERWIG ones. The only exception to this are the limits for the s generation, 
which are given default values by HERBVI as described in section 3.5. 



4 Conclusion 

We have described the various elements of the HERBVI program. This package takes 
the form of a hard scattering process event generator that interfaces to HERWIG, which 
is a general-purpose Monte Carlo program for simulating particle interactions involving 
hadrons. 

The phenomenological model which forms the motivation behind the HERBVI package 
is based upon the features of instanton-induced BLNV calculations, and we have briefly 
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HERBVI routine names 


Routine 


Function 


HVANAL 


Initial analysis 


HVCBVI 


Final state clustering 


HVCBVT 


Find hard parton parent 


HVETWO 


2 ^ 2 + n hard process 


HVHBVI 


HERWIG interface routine 


HVHEVT 


BLNV weight generation 


HVHGEN 


BLNV event data 


HVHINI 


Initialise derived parameters 


HVHMWE 


Multi-W weight generation 


HVHMWG 


Multi-W event data 


HVHPQM 


CKM matrix projection 


HVINIT 


Main initialisation 


HVRBIN 


Generate binomial distribution 


HVRBOS 


Generate boson distribution 


HVSFUN 


Generate parton x distribution 


HVSGEN 


Generate s distribution 



Table 3: Main HERBVI routines. 



HERBVI control flags 


Flag 


.TRUE. 


.FALSE. 


Comments 


HVCONT(l) 


Unweighted 


User weighted 


Phase space generation 


HVC0NT(2) 


BLNV 


Multi-W 


(internal flag) 


HVC0NT(3) 


On 


Off 


CKM matrix projection 


HVC0NT(4) 


On 


Off 


Uniform s generation 


HVC0NT(5) 


Off 


On 


Cut-off on s generation 


HVC0NT(6) 


Off 


On 


User-supplied s generation 


HVC0NT(7) 


On 


Off 


Set ub — 30 


HVC0NT(8) 


Off 


On 


User-supplied ub generation 


HVCONTO) 


On 


Off 


Allow Higgs production 



Table 4: Control flags. The default setting for each flag is . TRUE . , as set by HVINIT. 
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reviewed its features. Of particular importance from an event generation point of view is 
the parametrically large number of gauge bosons associated with a BLNV event, and we have 
addressed the issue of efficient phase space generation for such outgoing final states. The 
simulation of high-multiplicity gauge boson production with conservation of baryon- and 
lepton- number, which is expected to be the principal background to the BLNV processes 
considered here, is also implemented. 

The violation of baryon number requires particular attention when interfacing to a 
general simulation package such as HERWIG. Specifically, the formation of the extra 
(anti)baryons that violate baryon number conservation requires the correct colour con- 
nections between partons to be made after the emission of QCD radiation from the partons 
involved in the hard scattering. The interface between the HERBVI package and HERWIG 
allows these connections to be made at the appropriate time during event generation. 

Finally, we have described the provision that has been made within HERBVI to allow 
users to provide their own distributions. The distributions contained within the package 
correspond to the models discussed in 0. The interface to user-supplied distributions has 
been constructed in a straightforward manner to enable the study of BLNV processes with 
other parton-level cross sections and boson multiplicity distributions. 
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Appendix A: Installing HERBVI 

The HERBVI package is designed to be used with the HERWIG Monte Carlo simulation 
package. Versions 5.7 and later of HERWIG contain hooks that allow the two packages 
to link together in a relatively simple manner. There are two steps required to link the 
packages. 

• The dummy subroutines HVCBVI and HVHBVI must be deleted from the HERWIG 
package. Failure to do this will cause the program to terminate with an error message 
when the first attempt to generate a BLNV event is made. 

• The parameters NMXHEP, NMXPAR, MODMAX and NMXJET must be increased from their 
default values. The default values are set in the file HERWIGxy . INC, where x.y is the 
program version number, and the (recommended) new values are shown in Table |[ 
If this change is not made then HERWIG will run out of internal working space with 
unpredictable results. 

The first change can be made by deleting the dummy routines from the HERWIG 
source code. The second change must be made by editing the parameter settings in the 
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Parameter settings for HERBVI 


Parameter 


Default 


HERBVI 


NMXHEP 


2000 


15000 


NMXPAR 


500 


2000 


MODMAX 


5 


200 


NMXJET 


200 


500 



Table 5: Recommended HERWIG parameter settings for use with the HERBVI package. 

HERWIGxy . INC file. The changes in these parameters are required to give the HERWIG 
package enough working space to process the high-multiphcity BLNV events. Note that the 
formatting of printed output from HERWIG is automatically changed for NMXHEP values 
greater than 9999, in order to ensure that the 80-column output is readable. 

It should be stressed that, after the above changes have been made, the HERWIG 
package must be recompiled to reflect the change in size of the common blocks. 

Finally, we remind the reader that to use HERBVI with user-defined distributions 
(HVCONT (6) and/or HVCONT (8) = . FALSE . ) , the appropriate subroutines HUSGEN and/or HURBOS 
must be supplied by the user and the corresponding dummy subroutines deleted from the 
HERBVI package. 



Appendix B: Test run output 

The output of the HERBVI package, after processing by the HERWIG program, is contained 
as a list of final state hadrons and leptons within the standard HEPEVT common block. As 
there are typically ~ 10^~^ particles produced per HERBVI event, is is not practical to show 
one of these events here. Instead, we show the output of the test program supplied with 
the HERBVI package. This program generates 1000 events, with a total collider energy of 
40 TeV, using the IPR0C=7200 option. 

The implementation of random number generators is, generally speaking, not machine- 
independent due to numerical precision effects. Therefore, one should expect to get similar, 
but not necessarily identical, results when comparing the output of HERBVI runs from 
different computer architectures. 



OUTPUT ON ELEMENTARY PROCESS 



NUMBER OF EVENTS 
NUMBER OF WEIGHTS 
MEAN VALUE OF WGT 
RMS SPREAD IN WGT 
ACTUAL MAX WEIGHT 
ASSUMED MAX WEIGHT 



992 
52930 
1 . 6247E-04 
6 . 0487E-04 
7.7377E-03 
8 . 2604E-03 
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PROCESS CODE IPROC = 
CROSS SECTION (PB) = 
ERROR IN C-S (PB) = 
EFFICIENCY PERCENT = 



7200 

0.1625 
2.6291E-03 
1.967 



HERBVI parameters 



Total collider centre-of-mass energy (GeV) 
Minimum root s value (GeV) 
Maximum root s value (GeV) 
Parton-level cross section (nb) 
Number of gauge bosons set to fixed value 
Distribution of Higgs boson number 



0.40000E+05 

0.19596E+05 
0.35777E+05 
O.lOOOOE+01 
30 
Default 



Number of HVCBVI errors 
Number of events 

Percentage of events with BVI rej 
Number of BVI events 
Number of Multi-W events 



1000 
0.800 
502 
498 
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